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ABSTRACT 

An efficient search algorithm is described for the location of equilib- 
rium states in a search set of states which differ from one another only by 
the choice of pure phases. The algorithm has three important characteris- 
tics: (1) it ignores states which have little prospect for being an im- 

proved approximation to the true equilibrium state; (2) it avoids states 
which lead to singular iteration equations; (3) it furnishes a search his- 
tory which can provide clues to alternative search paths. 


INTRODUCTION 

The numerical computation of the composition and other properties of 
equilibrium states has been a continuing challenge for many years. As 
computers have become more powerful, people have attempted larger and more 
complex equilibrium calculations. These more complex problems naturally 
require greater amounts of computer time but, more importantly, they in- 
crease the probability that numerical obstacles to successful computations 
will arise. The potential impediments to success include inadequate initial 
estimates, convergence problems, difficulties with numerical significance and 
singular iteration equations. Difficulties such as these are often interre- 
lated and more than one of them may arise to plague the unwary thermodynami- 
cist in a particular problem. It is desirable, insofar as possible, to an- 
ticipate and avoid numerical problems to enhance the probability of success. 

Let me briefly review the nature of equilibrium computations. An equi- 
librium state corresponds to an extremum in a function over some selected 
class of candidate states; for example, it might be regarded as a maximum in 
the entropy or a minimum in a free energy. Thus, in principle, the location 
of an equilibrium state can be regarded as a two step procedure: the calcu- 

lation of the thermodynamic properties of all candidate states followed by a 
comparison of their entropies (or free energies). The greater the number of 
candidate states, the more difficult and time consuming the computational 
process and in short order it becomes quite impractical. The number of can- 
didate states usually grows with the number of chemical species being consi- 
dered. Hence an obvious way to limit the number of states is to restrict the 
number of species. But sometimes this may be undesirable because in so doing 
we might unwittingly eliminate the true equilibrium state from consideration. 


*Presented, in part, at Calphad XI held at the Argonne National Laboratory, 
Argonne, Illinois, May 16-20, 1982. 
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There is an alternative and generally more practical approach to the lo- 
cation of equilibrium states when there are many states to consider. One 
can begin at some initial state, explore for a state of higher entropy in 
the vicinity of this state and then step to an adjacent state of greater 
entropy. This permits one to ignore states of lower entropy and insures 
that each successful step will yield a state of higher entropy. The chief 
disadvantage with this method is that one cannot be certain that a partic- 
ular starting point will produce a search path that will ultimately reach 
the true equilibrium state; it could just as easily terminate at a local 
entropy maximum or terminate prematurely because of numerical difficulties. 
Further, even if the path were to terminate at the true equilibrium state 
one might not recognize it as the global maximum. A partial resolution of 
this dilemma is provided by trying different starting points and comparing 
the entropies at the end points of the various paths. This is illustrated 
schematically in fig. 1 for two search paths. In this illustration the 
first path gives a better approximation to the equilibrium state than the 
second but neither has located the true equilibrium state. 

Clearly the prospect of a successful search is enhanced by increasing 
the number of successful steps in the search path beginning at a given 
point. The number of steps can be increased by assiduously avoiding both 
states which are likely to create numerical problems and states which have 
little prospect for being the equilibrium state. I shall describe a search 
algorithm of this type for an important class of thermodynamic states; those 
obtained from a particular state by the addition of one or more pure phases 
in equilibrium. The algorithm has three beneficial characteristics: (1) it 

ignores states which have little prospect for yielding an improved approxi- 
mation to the true equilibrium state; (2) it avoids states which lead to 
singular iteration equations; (3) it furnishes a search history which might 
provide clues to~ alternative search paths. 


EQUILIBRIUM EQUATIONS 


Before I can describe the search algorithm it is necessary to look at 
the equation which must be satisfied by an equilibrium state. The governing 
equations are a set of nonlinear algebraic equations (Zeleznik and Gordon, 
1968; Zeleznik, 1981) obtained from a constrained extremum calculation. 


j = 1,2,..., m^ , ot — 1,2,. ..,p 
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i = 1 ,2, . . . 


In (1) the Greek index a is a phase index used to label each of the p 
phases, the Latin index j labels the ma species within a phase and the 
Latin index i labels the sl composition constraint relations. The quantity 


vi“ is the chemical potential of species j in phase a and is a function of 
the phase composition. The x-j are Lagrangian multipliers associated with 
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each of the 4 constraint relations. The doubly subscripted quantity 

is the set of stoichiometric coefficients representing the moles of component 
i (i^^ independent species) in one mole of the species in phase a. 

Finally n“ is the moles per unit mass of species j in phase a and gives 

the moles of component i per unit mass. Suppose that the phases Pq P 

are pure phases. Then for Pq 1 « <. P» = 1 cind u“ is just the chemTcal 

potential of the pure species u* which is independent of composition. Thus, 

in this case, ( 1 ) contains a subset of (1 + p - pq) linear equations which 
are independent of the composition variables. 


4 

^ a = Pq, Pq + 1,..., p 


( 2 ) 


Naturally any solution of (1) must provide a solution for (2). Conversely, 
if there is no solution of (2) there can be no solution for (1). Consequent- 
ly, it is important to focus on those states for which we can be relatively 
certain that ( 2 ) will possess solutions as a means for extending the search 
path. Fortunately, the properties of linear systems of equations are well 
understood and this understanding can be applied advantageously to the par- 
ticular set of equations shown in (2). The relevant ideas are most easily 
illustrated with a simple numerical example. 

A system of linear equations such as (2) can be written most compactly 
in matrix notation as 


Ax = c (3) 

where A is a rectangular matrix, c is a known column vector and x an 
unknown column vector. Linear equations are classified as either possessing 
solutions (consistent equations) or not possessing solutions (inconsistent 
equations). Consistent equations may possess one or many solutions. The 
matrix A and the vector c together determine whether the equations are 
consistent or inconsistent. For example let A be a 3 X 4 rectangular 
matrix. 

( 1 2 3 4\ 

2 3 4 6 ) 

3 5 7 loy 

Then x, if it exists, is a four component vector x' = (xj^, X 2 , X 3 , X 4 ) 
and c is a three component vector c' = (c^, 03 , C3) where a prime denotes 
the transpose of a matrix or a vector. Suppose that c' = 6 ' = ( 0 , 1, 1). 

Then it may be readily verified that x' = x' = (1, 1, -1, 0) is a solution 
of (3) and hence, for this choice of A and c, the equations are consistent. 
As a matter of fact there are infinitely many solutions. To demonstrate this 
we only need to verify that the two vectors 
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satisfy Ax? = 0 = Ax^. Such vectors are said to be in the null space of A. 
i. 0 0 

Thus the vector x = x + b^x^ + is a solution for arbitrary numbers 

and B 2 * Actually this is the general solution to the problem. On the 

other hand suppose c' = c' s (1, 1, 1). Then it turns that (3) becomes an 
inconsistent set of equations because if a solution existed it would lead to 
an inconsistency. To see that this is true we observe that if w'h(1, 1, -1) 
then w'A = 0. But if a solution exists then 

1 = (1 + 1 - 1) = w‘c = w'Ax = 0‘x = 0 

which is clearly inconsistent. A vector such as w is said to be in the left 
null space of A or, equivalently, in the null space of A*. Thus we see 
that (3) is solvable for x if and only if c satisfies the condition 

w'c = C]^ + C 2 - C 3 = 0 

In general. Ax = c is solvable if and only if w'c = 0 for each indepen- 
dent vector w in the left null space of A. 

The ideas which were illustrated by example in the preceding paragraph 
can be established quite generally and are done so in texts on linear alge- 
bra. For our purposes the important feature is that each independent vector 
in the left null space of A imposes an independent constraint which must 
be satisfied by the vector c if Ax = c is to be consistent. Each addi- 
tional constraint reduces the likelihood that a given vector c will permit 
a solution. However, if there are no vectors in the left null space of A 
then a solution will always exist. If the relation w'A =0 is written out 
in terms of the components then it can be recognized as implying that the 
rows of A are linearly dependent. Thus if the rows of A are linearly 
independent then there will be no vector in the left null space of A and 
no constraints imposed on c. 

Suppose we now choose A to be the coefficient matrix of (2). Then its 
rows are 


/ a a a \ j. 1 

’ '*21 ’"**'*)ll a — Pq»Pq !»•••» P 

and each row is composed of the stoichiometric coefficients for the species 
which forms one of the 1 + p - pg pure phases. So long as the pure 
phase species are independent of each other there will be no restrictions 
imposed on the permissible values for the right side of (2). Thus, from the 
previous discussion, it seems desirable to confine a search path to states 
whose pure phases are independent species. This does not mean that depen- 
dent pure phases cannot coexist. For example, on any typical winter day one 
can find liquid water and ice coexisting and these are certainly dependent 
species since they have identical chemical formulas. It only means that a 
state containing both water phases is less likely to arise as an equilibrium 
state than one containing only one of them. The initial phases of a search 
for an equilibrium state should be confined to the more likely states. Those 
states containing dependent pure phases can be searched after all of the in- 
dependent pure phase states have been explored. As a practical matter candi- 
date states with dependent pure phases may have to be treated as special 
cases rather than as part of a general search procedure. 
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We can now incorporate the ideas from the foregoing discussion in a gen- 
eral search algorithm over a set of thermodynamic states which differ from 
one another only in the choice of independent pure phases. 


Step 1. 
Step 2. 
Step 3. 
Step 4. 

Step 5. 


Step 6. 


Step 7. 


Specify an initial state with p independent pure phases. 
Specify a search set of pure species. 

Calculate properties of the state. 

If any pure phases are negative remove all of them and go to 
Step 3. 

For each species in the search set calculate 

i 






where p. is the chemical potential of the 

J 


J 

i=l 

jth species and vjj is its set of stoichiometric coefficients. 
Select as the (p+lj^^ pure phase a species from the search set for 
which ^ < 0, ^ is a minimum over the search set and such that the 
species is independent of the p phases already present. 

If a (p+1)^” pure phase is located go to Step 3, otherwise terminate 
the search. 


The criterion 0 < 0 for pure phase inclusion. Step 6, is used in the current 
NASA equilibrium program and was also used by its ancestral programs (Gordon 
and McBride, 1976; Gordon and Zeleznik, 1963; Zeleznik and Gordon, 1962). The 
simultaneous removal, in Step 4, of all pure phases with negative moles differs 
from the procedure used in the NASA equilibrium programs where only one such 
phase is removed. The reason for this change is that it introduces an asymme- 
try in the search path which should help to minimize the occurrence of closed 
loop paths. Occasionally it may also accelerate progress toward the final 
state. The requirement for independent species in Step 6 has a twofold pur- 
pose. First, it restricts the search to independent pure phases which have 
the greatest prospect for yielding a solution. Second, it avoids the class of 
singular iteration equations associated with dependent pure phases. 


AN APPLICATION OF THE SEARCH ALGORITHM 

I shall illustrate the use of the search algorithm by applying it to a 
large, difficult problem that arose in a study of the corrosion of turbine 
blade materials by deposits from high temperature combustion gases (Santoro 
et al., 1981, 1980). The problem involves 19 chemical elements present in 
widely varying amounts as is shown in Table 1. Further, because many of the 
species of interest are present in trace amounts it was desirable to calculate 
accurately mole fractions as small as lO'^O, jhe large disparity in elemen- 
tal amounts and the requirement for high accuracy in the mole fractions pre- 
sented a severe challenge to the numerics. Even with double precision arith- 
metic (approximately 14 digits) the results are just above the numerical noise. 
Calculations were carried out at one atmosphere pressure and a temperature of 
1000 K. The initial states were always composed of the gaseous phase plus var- 
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ious pure phases. Nearly 500 different chemical species were considered in 
the equilibrium computations. 

Table 2 shows the initial and final states for six different search paths, 
the number of steps in each path and the Gibbs free energy of the final state. 

Paths 3 and 4 terminated at the same state while paths 5 and 6 also terminated 
at identical states. The common terminus of paths 3 and 4 represents the best 
estimate of the equilibrium state as judged by the value of the Gibbs free en- 
ergy. The paths themselves are schematically presented in Table 3 where "+" 
denotes an addition of a pure phase, denotes a deletion of a pure phase 
and "±" represents a null step. A null step means that the algorithm would 
have added a pure phase but could not do so because it would have been linear- 
ly dependent on the pure phases already present. Paths 3 and 5 are not shown. 

Path 3 has no steps and path 5 has two steps, the removal of Al2Si05 fol- 
lowed by the addition of CaAl407. 

Path 1 begins at an initial state containing only the gaseous solution and 
terminates in a tight three step loop without encountering any null steps. 

Table 4 shows the properties of states encountered by path 1 . The initial 
state for path 2 is a modification of the final state for path 1 . Neither 
Al2Ti05 nor CaaTi^Oj^Q were chosen for the initial state of path 2 because they 
were involved in the tight loop terminating path 1 . The species Ti^OT was re- 
placed by Ti02 and AlgSi20i3 was replaced by AI0O3 and Si02 with whicn it is 
linearly dependent. Finally, Fe304 was replaced By Fe203 and the species 
CaS04, Na2S04, PbS04 and ZnO were added. The ensuing search path contains two 
null steps; one involving Al2Si0c to terminate the search and another involving 
Al2Ti05 at step 2 . Although Al2Ti05 was involved in a null step, it was 
subsequently added to the system at step 8. 

The initial states for paths 3 and 4 are derivatives of the final state for 
path 2 . Since path 2 was completed with an Al2Si0q null step it suggests that new 
paths be tried which include it as a pure phase. Thus, for path 3 , AlgSi2033 was 
omitted from consideration in the search and was replaced by Al2Si0g in the initial 
state. This search terminated with no additional steps being taken. For path 4 , 
AlgSi2033 was again omitted from consideration and was replaced in the initial state 
by the combination of AI2O3 and Si02. 

Like paths 3 and 4 , the initial states for 5 and 6 are closely related to 
the final state for 2 . Rather than omitting Al6Si20l3, it was AI2O3 which 
was excluded from the four species dependent set of AI2O3, SiO^, AloSiOg 
and AlgSi20i3. For path 5 the species AlgSio033 was supplanted by Al2Si0g 
in the initial state. For path 6 it was replaced by Si02- 


CONCLUSIONS 

An algorithm has been described which offers considerable assistance in 
the determination of equilibrium states containing pure phases. Its utility 
increases with the complexity of the system being calculated. 
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TABLE 1. 



- DISTRIBUTION OF ELEMENTS IN THE SAMPLE PROBLEM. 


c 

KG-ATOMS/KG 
0.2796061 lE-02 

H 

0.40735049E-02 

N 

0.51352555E-01 

0 

0.14417322E-01 

s 

0.35989835E-05 

Ar 

0.30643626E-03 

V 

0.37750290E-09 

Ca 

0.95960942E-08 

Mg 

0.15824376E-07 

A1 

0.71273032E-07 

Cl 

0.54242432E-07 

Na 

0.83648392E-08 

K 

0.49 185 178E-08 

Si 

0.6847151 IE-08 

Fe 

0. 3443437 7E -08 

2n 

0.29413516E-08 

Pb 

0.9281161 IE-09 

Ba 

0.14003181E-08 

Ti 

0.40147322E-08 



SEARCH PATH 


QE POOR QUALITY 


TABLE 2. - PURE PHASE SPECIES IN INITIAL AND FINAL STATES FOR SEARCH PATHS 


PURE PHASES INITIAL FINAL INITIAL FINAL INITIAL FINAL INITIAL FINAL INITIAL FINAL INITIAL FINAL 

AI 2 O 3 XXXXXXO 0 

SiO^ XXX 

Al 2 Si 05 XX XX 

^ X 0 0 XX X 


GIBBS ENERGY -2156.7272 6236 -2156.7275 3378 
PER UNIT MASS' 


-2156.7275 3734 


SEARCH 


NOTE: X PURE PHASE PRESENT 

0 PURE PHASE OMITTED FROM CONSIDERATION IN THE 


•2156.7274 7210 


9 




TABLE 3. - SEARCH PATH DETAILS 
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TABLE 4. - STATES ALONG SEARCH PATH 1 

THERnOOYNAniC EQUUXBRIUH PROPERTIES AT ASSIGNED 
TEMPERATURE AND PRESSURE 



CHEMICAL 1 

FORMULA 




(SEE NOTE) 

CAL/MOL 


DEG K 

C/CC 

FUEL 

C 

\ 

00000 

H 

1 

07090 N 

0.00890 

0 0.03360 

0 9968<^0 

631 

000 

L 

298 

15 

0 

0000 

fUU 

S 

1 

00000 






0 003000 

-17 

309 

5 

298 

15 

0 

0000 

OXIDANT 

M 

1 . 

S6I76 

0 

0 

A1959 AR 

0.00932 

C 0.00030 

0 990000 

-28 

200 

G 

298 

15 

0 

0 000 

OXIDANT 

H 

2 . 

00000 

0 

1 

00000 



0 010000 

-57795 

926 

G 

298 

15 

0 

0000 

ruci 

V 

1 

00000 






0000005 

0. 

000 

S 

298 

15 

0 

0000 

FULL 

CA 

1 . 

00000 






OOOOIDO 

0 

000 

S 

298 

15 

0 

0000 

f UFL 

MG 

1 

00000 






0000100 

0 

000 

S 

298 

15 

0 

0000 

FUEL 

AL 

1 

00000 






.0000500 

0 

000 

S 

298 

15 

0 

0000 

1 UFL 

Cl 

2 

00000 






.0000500 

0 

000 

G 

298 

15 

0 

0000 

FUFL 

NA 

1 

00000 






.0000050 

0 

000 

S 

298 

15 

0 

0000 

FUF.L 

K 

1 

00000 






0000050 

0 

000 

s 

298 

15 

0 

0000 

FUF L 

SI 

1 

00000 






.0000050 

0 

000 

s 

298 

15 

0 

0000 

FUEL 

FE 

1 

00000 






0000050 

58 

876 

s 

298 

15 

0 

0000 

FUEL 

ZN 

1 

00000 






.0000050 

0 

000 

s 

298 

15 

0 

0000 

fuel 

PB 

1 

00000 






0000050 

0 

000 

s 

298 

15 

0 

0000 

FUFL 

BA 

1 

00000 






OOOOOSD 

0 

000 

s 

298 

15 

0 

0000 

FUEL 

TI 

1 

00000 






.0000050 

0 

000 

s 

298 

15 

0 

0000 

0/F = 

25 0000 

PERCENT FUELS 

3 0^62 

EQUIVALENCE RATIOS 0 5297 

PHIS 0 5076 

REACTANT 

DENSITYs 

0 . 

0000 


IMERMODYHAMIC PROPERTIES 


p , 

ATM 

1 0000 

1 0000 

1 0000 

1 .0000 

1 0000 

1 0000 

1 0000 

1 0000 

1 0000 

1 0000 

1 0000 

1 0000 

1 0000 

] . 

Drr. K 

1000 0 

1000 0 

1000 0 

1000 0 

1000.0 

1000 0 

1000 0 

1000 0 

1000 0 

lOOO 0 

1000 0 

1 000 0 

1000 0 

RMO. G/CC 

3 5629-*^ 

3 5629-9 

3 5629-9 

3.5629-9 

3 5629-9 

3.5629-9 

3 5629-9 

5 5629-9 

3 5629-9 

5 5629-9 

3 5629-9 

3 5629-9 

3 5629-9 

H . 

CAl/G 

-199 77 

-199 77 

-199 77 

-199 77 

-199 78 

-199 78 

-199 78 

-199.78 

-199 78 

-199 78 

-199 78 

-199 78 

-199 78 

u. 

CAl /G 

-262 75 

-262 75 

-262 75 

-262 75 

-262.76 

-262 76 

-262 76 

-262 76 

-262 76 

-262 76 

-262 76 

-262 76 

-262 76 

G. 

CAL/G 

-2156 72 

-2156 72 

-2156 72 

-2156 . 72 

-2156 73 

-2156 73 

-2156 73 

-2156 73 

-2156 73 

-2156 73 

-2156 73 

-2156 73 

-2156 73 

b . 

CAL/(G)(K) 

1 9620 

1 9620 

1 9620 

1 9619 

1 9619 

1 9619 

1.9619 

1 9619 

1 . 9619 

1 9619 

1 9619 

1 9619 

1 9619 

M. 

not UT 

29 232 

29 232 

29 232 

29 232 

29 252 

29 232 

29 252 

29.232 

29.232 

29 232 

29 232 

29 232 

29 232 

(DlV/DlPiT 

-1 00001 

-1 00001 

-1 OOOCl 

-1 00001 

-1 00001 

-1 00001 

-1 00001 

-1 00001 

-1 00001 

-1 00001 

-1 00001 

-1 OOOOl 

-1 00001 

(DLV/DLT)P 

1 0001 

1 00 0 1 

1 0001 

1 0 001 

1 0001 

1 0001 

1 0001 

1 0001 

1 0001 

1 000 1 

1 00 0 1 

1 0001 

1 0001 

Cl' . 

CAl/IGHK) 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

0 2873 

GAMMA (S) 

1 3101 

1 3101 

1 3101 

1 3101 

1 3101 

1 3101 

1.3101 

1 3101 

1.3101 

1 3101 

1 3101 

1 3101 

1 3101 


MUIE FRACTIONS 
SOLUTION 1 


CH202 

1 

090-18 

1 

090-18 

1 

090-18 

1 

090-18 

1 

090-18 

1 

090-18 

1 090-18 

1 

090-18 

2 

.090-28 

1 

090-18 

1 

.090-18 

1 

090-18 

1 

090-18 

K2C03 

8 

863-19 

8 

863-19 

8 

803-19 

8 

769-19 

a 

769-19 

8 

757-19 

8.803-19 

8 

753-19 

8 

752-19 

8 

752-19 

8 

706-19 

8 

758-19 

8 

706-19 

ZM 

1 

919-12 

1 

919-12 

1 

.919-12 

1 

919-12 

1 

919-12 

1 

919-12 

1.919-12 

1 

919-12 

1 

919-12 

1 

919-12 

1 

.919-12 

I 

919-12 

1 

.919-12 

ZMO 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8.5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

8 

5980-8 

Air 1 

9 

598-20 

9 

296-20 

9 

397-20 

1 

867-20 

8 

912-21 

8 

969-21 

9 987-36 

9 

289-33 

9 

291-35 

9 

291-35 

9 

911-35 

9 

618-36 

9 

911-55 

A 1 C L 2 

J 

912-16 

1 

766-16 

1 

899-16 

8 

199-17 

3 

.699-17 

3 

.795-17 

1 933-32 

1 

900-29 

1 

902-29 

1 

905-29 

2 

010-29 

2 

093-32 

2 

010-29 

AICL3 

1 

067-10 

9 

863-11 

1 

095-10 

9 

833-11 

2 

178-11 

2 

223-11 

1 118-26 

I 

150-23 

1 

152-23 

1 

152-23 

1 

.229-23 

1 

219-26 

1 

229-23 

AlOCL 

3 

117 -9 

2 

879 -9 

2 

981 -9 

1 

266 -9 

5.703-10 

5 

.792-10 

3 092-25 

2 

908-22 

2 

909-22 

2 

909-22 

2 

990-22 

3 

131-25 

2 

990-22 

ALOH 

2 

551-17 

2 

355-17 

2 

359-17 

9 

808-13 

9 

918-18 

9 

918-18 

2.902-33 

2 

253-50 

2 

255-30 

2 

235-50 

2 

233-30 

2 

908-35 

2 

233-30 

AL02 

9 

855-12 

9 

982-12 

9 

982-12 

1 

867-12 

8 

910-13 

8 

910-13 

9.573-28 

9 

250-25 

9 

250-25 

9 

250-25 

9 

250-25 

9 

583-28 

9 

250-25 

Al 02H 

2 

080 -6 

1 

920 -6 

1 

920 -6 

8 

.000 -7 


603 -7 

3 

603 -7 

1 959-22 

1 

821-19 

1 

821-19 

1 

.821-19 

1 

821-19 

1 

969-22 

1 

821-19 

AL202 

7 

072-20 

6 

027-20 

6 

026-20 

1 

. 096-20 

2 

122-21 

2 

122-21 

0 000 0 

0 

000 0 

0 

000 0 

0 

000 0 

1 

*>00 0 

0 

000 0 

0 

000 0 

AR 

8 

95/7-3 

8 

9577-3 

8 

9577-3 

8 

9578-3 

8 

9578-3 

8 

9578-3 

8 9578-3 

8 

9578-3 

8 

9578-5 

8 

9578-3 

8 

9578-3 

8 

9578-3 

8 

9578-3 

B*.CL 

9 

995-19 

9 

995-19 

9 

991-19 

9 

936-19 


936-19 

«. 

939-. 9 

9 991-19 

9 

939-19 

9 

939-19 

9 

955-19 

1 

.587-31 

1 

591-31 

1 

587-31 

DACL2 

2 

053 -8 

2 

059 -8 

2 

126 -8 

2 

165 -8 

~ 

165 -5 

2 

.179 -8 

2.126 -8 

2 

183 -8 

2 

189 -8 

2 

189 -8 

7 

222-21 

6 

809-21 

7 

222-21 

6AOH 

5 

209-18 

5 

207-18 

5 

029-18 

9 

925-18 

9 

929-18 

9 

869-18 

5.029-18 

9 

876-18 

9 

876-18 

9 

879-18 

1 

526-30 

1 

526-50 

1 

526-30 

n/‘ 0 .‘*M 2 

2 

090 -8 

2 

039 -8 

1 

968 -8 

1 

929 -8 

1 

928 -8 

1 

.915 -8 

1 968 -8 

1 

910 -8 

1 

.909 -8 

1 

909 -8 

5 

975-21 

5 

975-21 

5 

975-21 

CO 

1 

562-U 

1 

562-11 

1 

562-U 

1 

562-11 

1 

562-U 

1 

562-U 

1 562-11 

1 

562-11 

1 

.562-11 

1 

562-U 

1 

562-U 

1 

562-U 

1 

562-U 

C02 

8 

1739-2 

8 

1739-2 

8 

1735-2 

8 

1735-2 

8 

1735-2 

8 

1735-2 

8 1735-2 

8 

.1735-2 

8 

1755-2 

8 

1735-2 

8 

1735-2 

8 

1735-2 

8 

1735-2 

CACL 

2 

095-18 

2 

096-18 

9 

268-19 

0 

000 0 

0 

000 0 

0 

. 000 0 

1 059-18 

1 

953-30 

1 

959-50 

1 

959-30 

1 

999-30 

1 

073-18 

1 

999-30 

CACL2 

9 

198 -8 

9 

201 -8 

1 

971 -8 

0 

0 00 0 

0 

000 0 

0 

000 0 

2 291 -6 

3 

178-20 

3 

181-20 

3 

182-20 

3 

361-20 

2 

399 -8 

3 

361-20 

tAOM 

3 

678-17 

3 

678-17 

1 

609-17 

0 

000 0 

0 

000 0 

0 

000 0 

1 829-17 

2 

952-29 

2 

952-29 

2 

952-29 

2. 

,952-29 

1 

819-17 

2 

952-29 

CAnZH2 

2 

385 -7 

2 

385 -7 

1 

093 -7 

0 

000 0 

0 

000 0 

0 

000 0 

1.186 -7 

1 

590-19 

1 

590-19 

1 

590-19 

1 

590-19 

1 

176 -7 

1 

590-19 

CL 

6 

986-10 

6 

909-10 

7 

235-10 

7 

378-10 

7 

. 379-10 

7 

929-10 

7 238-10 

7 

995-10 

7 

998-10 

7 

999-10 

7 

656-10 

7 

939-10 

7 

656-10 

Cl n 

1 

379-n 

1 

375-11 

1 

929-11 

1 

951-11 

1 

952-11 

1 

962-11 

1.929-11 

1 

965-11 

1 

965-11 

1 

966-11 

1 

506-11 

1 

963-11 

1 

506-11 

Cl 0.’ 

9 

570-16 

9 

572-16 

9 

755-16 


826-16 

9 

827-16 

A 

860-16 

9.735-16 

9 

870-16 

9 

872-16 

9 

873-16 

5 

008-16 

9 

863-16 

5 

008-16 

U 2 

3 

259-12 

3 

262-12 

3 

999-12 

3 

635-12 

3 

.636-12 

3 

686-12 

5 999-12 

3 

701-12 

3 

709-12 

3 

706-12 

3 

915-12 

3 

691-12 

3 

915-12 

CL20 

2 

779-20 

2 

781-20 

2 

983-20 

3 

099-20 

3 

100-20 

3 

192-20 

2 985-20 

3 

155-20 

3 

158-20 

3 

159-20 

3 

337-20 

3 

196-20 

3 

337-20 

rrci 2 

1 

977-10 

1 

978-10 

1 

585-10 

1 

697-10 

1 

697-10 

1 

.670-10 

1 585-10 

1 

677-10 

1 

678-10 

1 

626-15 

1 

718-15 

1 

620-15 

1 

718-15 

ri ct 3 

7 

599-19 

7 

559-19 

6 

389-19 

8 

883-19 

8 

888-19 

9 

070-19 

8 389-19 

9 

127-19 

9 

138-19 

8 

857-19 

9 

616-19 

8 

803-19 

9 

616-19 

n 0 

6 

999-17 

6 

999-17 

6 

999-17 

6 

998-17 

6 

998-17 

6 

998-17 

6 999-17 

6 

998-17 

6 

998-17 

6 

299-22 

6 

299-22 

6 

299-22 

6 

299-22 

FFO'»H2 

1 

005 -7 

1 

005 -7 

1 

005 -7 

1 

005 -7 

1 

.005 -7 

1 

005 -7 

1.005 -7 

1 

005 -7 

1 

005 -7 

9 

759-13 

9 

739-13 

9, 

,739-13 

9 

739-13 

FI.2CL9 

2 

519-19 

2 

523-19 

2 

901-19 

3 

131-19 

3 

133-19 

3 

219-19 

2 901-19 

3 

296-19 

3 

251-19 

3 

053-29 

3 

907-29 

3 

028-29 

3 

907-29 

M 

9 

129-15 

9 

129-15 

9 

125-15 

9 

125-15 

9 

125-15 

9 

125-15 

9.125-15 

9 

125-15 

9 

125-15 

9 

125-15 

9 

125-15 

9, 

,125-15 

9 

125-15 

MFL 

1 

3956-6 

1 

3961-6 

1 

3991-6 

1 

9210-6 

1 

9213-6 

1 

9309-6 

1.3991-6 

1 . 

.9539-6 

1 

9395-6 

1 

9398-6 

1 

9797-6 

1 

9319-6 

1 

9797-6 

UNO 

2 

699-19 

2 

699-19 

2 

699-19 

2 

699-19 

2 

699-19 

2 

699-19 

2.699-19 

2. 

,699-19 

2. 

699-19 

2. 

699-19 

2. 

,699-19 

2 

699-19 

2 

699-19 

HMD2 

1 

9519-9 

1 

9519-9 

1 

9519-9 

1 

9519-9 

1 

9519-9 

1. 

. 9519-9 

1,9519-9 

1 , 

.9519-9 

1 . 

9519-9 

1 

9519-9 

1 

9519-9 

1 . 

9519-9 

1 

9519-9 

Mi;03 

1 

906-11 

1 

906-11 

1 

906-11 

1 

906-11 

1 

906-11 

I. 

, 906-n 

1 906-11 

1 

906-11 

1 

906-11 

1 

906-11 

1 . 

.906-11 

1 

906-11 

1 . 

.906-11 

un » 

1 

9533-9 

1 

9539-9 

1 

9539-9 

1 

9539-9 

1 . 

,9539-9 

1 

9539-9 

1 9559-9 

1 

9539-9 

1 . 

9539-9 

1 . 

,9539-9 

1 

9539-9 

1 . 

.9539-9 

1 

9539-9 

M3 

1 

638-11 

1 

638-11 

1 

638-11 

1 

638-11 

1 

638-11 

1 

638-11 

I 638-11 

1. 

.638-11 

1 

638-11 

1 

638-11 

1 

638-11 

1 

638-U 

1 

638-11 

H20 

5 

9536-2 

5 

9536-2 

5 

9536-2 

5 

9537-2 

5 

9557-2 

5 

9537-2 

5 9537-2 

5 

9538-2 

5 

9558-2 

5 

9538-2 

5 

9538-2 

5 

9538-2 

5 

9538-2 

M202 

5 

559-11 

5 

559-11 

5 

559-11 

5 

559-11 

5 

559-11 

5 

559-11 

5.559-11 

5 

559-11 

5 

559-n 

5 

559-11 

5 

559-11 

5 

559-11 

5 

559-11 

M2S09 

9 

8726-9 

9 

8726-9 

9 

8727-9 

9 

8729-9 

9. 

8729-9 

9. 

8729-9 

9 8729-9 

9. 

8730-9 

9. 

8735-9 

9 

8735-9 

9 

8697-9 

9 

8697-9 

9 

8697-9 

K 

1 

090-16 

1 

039-16 

1 

036-16 

1 

039-16 

1 

039-16 

1 

035-16 

1.036-16 

1 

033-16 

1 . 

033-16 

1 

033-16 

1 

030-16 

1 

033-16 

1 

030-16 

KCl 

2 

9967-8 

2 

9976-8 

2 

5778-8 

2 

6225-8 

2 

6230-8 

2 

6390-8 

2.5779-8 

2 

6990-8 

2. 

6999-8 

2 

6953-8 

2 

7117-8 

2. 

6910-8 

2 

7U7-8 

KO 

1 

511-17 

1 

511-17 

1 . 

506-17 

1 

503-17 

1. 

503-17 

1 

502-17 

1,506-17 

1 . 

502-17 

1 

502-17 

1. 

502-17 

1 

998-17 

1 . 

502-17 

1 . 

998-17 

KOM 

1 

909-10 

1 

909-10 

1 

909-10 

1 

901-10 

1 . 

901-10 

1 . 

900-10 

1.909-10 

1 . 

900-10 

1 

900-10 

1 

900-10 

1. 

396-10 

1 

901-10 

1 . 

396-10 

K3CL2 

1 

958-12 

1 

959-12 

2 

087-12 

2 

160-12 

2 

161-12 

2 

187-12 

2 087-12 

2 

195-12 

2. 

197-12 

2 

198-12 

2 

309-12 

2 

190-12 

2 

309-12 

K/02H2 

2 

298-18 

2. 

298-18 

2 

282-18 

2 

279-18 

2 

275-18 

2. 

270-18 

2.282-18 

2 

269-18 

2. 

269-18 

2 

269-18 

2 

257-18 

2. 

271-18 

2 

257-18 

K.'»S09 

5 

9333-8 

5 

9329-8 

5 

8928-8 

5 

8709-8 

5. 

8702-8 

5 

8622-6 

5.8928-8 

5 

8597-8 

5. 

8593-8 

5 

8590-8 

5 

8259-8 

5 

8612-8 

5.8259-8 

Meet 

1 

007-18 

8 

327-19 

8 

621-19 

8 

787-19 

1 

063-18 

1 . 

070-18 

8.971-29 

1 

009-29 

1 

009-29 

1 

009-29 

1 

037-29 

8 

663-29 

1 

037-29 

nr,CL2 

1 

222 -9 

1 

012 -9 

1 . 

085 -9 

1 

127 -9 

1. 

363 -9 

1 . 

382 -9 

1 .066-19 

1 

305-20 

1 . 

306-20 

1 

307-20 

1 

381-20 

1 

119-19 

1 

381-20 

Mun 

2 

968-19 

2 

955-19 

2 . 

959-19 

2 

959-19 

2 

967-19 

2 . 

967-19 

2.9U-29 

2 

792-30 

2 

792-30 

2 . 

792-30 

2 

792-50 

2 

901-29 

2 

792-30 
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TABLE 4. - Concluded. 


ncoM 

7 

35S-16 

6 

082-16 

6 

081-16 

6.080-16 

7 353-16 

7.353-16 

5.975-21 

6 

917-27 

6 

.917-27 

6 

.917-27 

6 

917-27 

5 

949-21 

6 

.917-27 

MG0?H2 

A 

.61A -7 

3 

815 -7 

5 

81A -7 

3 81A -7 

A 612 -7 

A. 612 -7 

3 748-12 

4 

339-18 

A 

.339-18 

A 

.339-18 

4 

.339-18 

3 

732-12 

4 

-339-18 

M<3 

3 

39A'20 

3 

398-20 

3 

. 398-20 

3 398-20 

3 398-20 

3.398-20 

3 398-20 

3 

398-20 

3 

398-20 

3 

398-20 

3 

398-20 

3 

.398-20 

3 

.398-20 

MO 

2 

3612-S 

2 

3612-5 

2 

3612-5 

2 3612-5 

2.3612-5 

2.3612-5 

2 3612-5 

2 

3612-5 

2 

3612-5 

2 

3612-5 

2 

3612-5 

2 

3612-5 

2 

3612-5 

NOCL 

5 

03A-12 

3 

035-12 

5 

1A3-12 

3 203-12 

3.20A-12 

3 226-12 

3 143-12 

3 

233-12 

3 

.234-12 

3 

235-12 

3 

324-12 

3 

.228-12 

3 

324-12 

N02 

8 

. 57A8-7 

8 

57A9-7 

8 

.57A9-7 

8.57A9-7 

8.57A9-7 

8 57A9-7 

8.5749-7 

8 

5749-7 

8 

5749-7 

8 

5749-7 

8 

5749-7 

8 

5749-7 

8 

.5749-7 

H02CL 

1 

606-lS 

1 

607-15 

1 

66A-15 

1 696-15 

1.696-15 

1.708-15 

1 .664-15 

1 

711-15 

1 

712-15 

1 

713-15 

1 

760-15 

1 

709-15 

1 

760-15 

NOS 

6 

291-lA 

6 

.291-lA 

6 

.291-lA 

6. 291-lA 

6 291-lA 

6.291-14 

6.291-14 

6 

291-14 

6 

.291-14 

6 

.291-14 

6 

291-14 

6 

.291-14 

6 

.291-14 

N? 

7 

5056-1 

7 

5056-1 

7 

.5056-1 

7 5056-1 

7.5056-1 

7 5056-1 

7.5056-1 

7 

5056-1 

7 

5056-1 

7 

5056-1 

7 

5056-1 

7 

.5056-1 

7 

5056-1 

N20 

1 

5903-9 

1 

.5903-9 

1 

5903-9 

1.5903-9 

1.5903-9 

1.5903-9 

1.5903-9 

1 

5903-9 

1 

5903-9 

1 

5903-9 

1 

5903-9 

1 

.5903-9 

1 

.5903-9 

H20A 

8 

378-19 

8 

378-19 

8 

378-19 

8.378-19 

8.379-19 

8.379-19 

8.379-19 

8 

.379-19 

8 

.379-19 

8 

.379-19 

8 

379-19 

8 

379-19 

8 

579-19 

MA 

A 

696-16 

A 

696-16 

8 

679-16 

8.670-16 

a 670-16 

8.666-16 

8 679-16 

8 

665-16 

8 

.456-16 

8 

455-16 

8 

443-16 

8 

.458-16 

8 

443-16 

NACL 

2 

5299-8 

2 

.5309-8 

2 

6159-8 

2.663A-8 

2 6639-8 

2.6809-8 

2 6160-8 

2 

6862-8 

2 

6224-8 

2 

.6228-8 

2 

.6917-8 

2 

6183-8 

2 

6917-8 

NAO 

I 

930-16 

1 

930-16 

1 

927-16 

1.92A-16 

1 92A-16 

1 924-16 

1 927-16 

1 

923-16 

1 

.877-16 

1 

.877-16 

1 

.874-16 

1 

877-16 

1 

874-16 

NAOH 

1 

.A12-10 

1 

A12-10 

1 

AlO-lO 

1.A08-10 

1.A08-10 

1 408-10 

1.410-10 

1 

407-10 

1 

.373-10 

1 

.373-10 

1 

.371-10 

1 

374-10 

1 

.371-10 

NA2CL2 

A 

lAO-12 

A 

. 1A3-12 

A 

A26-12 

A 588-12 

A. 590-12 

4 649-12 

A 426-12 

A 

.667-12 

A 

448-12 

A 

449-12 

A 

.686-12 

A 

.434-12 

4 

686-12 

NA?02H2 

I 

8AA-17 

1 

.8AA-17 

1. 

.836-17 

1.832-17 

1.832-17 

1 831-17 

1.836-17 

1 

831-17 

1 

743-17 

1 

743-17 

1 

.738-17 

1 

744-17 

1 

738-17 

NA2SOA 

1 

095A-7 

1 

.0953-7 

1 , 

.0911-7 

1.0887-7 

1.0887-7 

1.0878-7 

1.0911-7 

1 

0875-7 

1. 

.0557-7 

1 

0356-7 

1 

.0321-7 

1 , 

.0358-7 

1 

0321-7 

0 

A 

916-11 

A 

S’6-11 

A 

916-11 

A.916-H 

A 916-11 

4.916-11 

4.916-11 

A 

916-11 

A 

916-11 

A 

.916-11 

A 

916-11 

4 

916-11 

4 

.916-11 

OH 

7 

6A50-8 

7 

eA50-8 

7 

6A50-8 

7.6A51-8 

7.6A51-8 

7 6451-8 

7.6451-8 

7 

6451-8 

7, 

.6451-8 

7, 

.6451-8 

7 

6451-8 

7 

6451-8 

7 

.6451-8 

02 

9 

.9081-2 

9 

.9081-2 

9 

.9081-2 

9.9081-2 

9.9081-2 

9.9081-2 

9.9081-2 

9 

9081-2 

9 

9081-2 

9 

.9081-2 

9 

.9081-2 

9, 

.9081-2 

9 

9081-2 

03 

2 

729-13 

2 

729-13 

2. 

.729-13 

2.729-13 

2.729-13 

2 729-13 

2.729-13 

2 

729-13 

2 

729-15 

2 

729-13 

2 

729-13 

2 

729-13 

2 

729-13 

> 3 

1 

200-12 

1 

.200-12 

1 

125-12 

1 085-12 

1.085-12 

1.070-12 

1.125-12 

1. 

.066-12 
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